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In this paper we use collocation methods for detecting blow-up solutions of nonlinear 
homogeneous Volterra-Hammerstein integral equations. To do this, first we generalize 
the results on existence and uniqueness given in a previous paper in order to consider also 
non locally bounded kernels, which are present in many cases of blow-up problems. Next, 
we introduce the concept of "blow-up collocation solution" and analyze numerically 
some blow-up time estimates using collocation methods in particular examples where 
the previous results about existence and uniqueness can be applied. Finally, we discuss 
the relationships between necessary conditions for blow-up of collocation solutions and 
exact solutions. 

1 Introduction 

Some engineering and industrial problems are described by explosive phenomena which are 
modeled by nonlinear integral equations whose solutions exhibit blow-up at finite time (see 
[H [2] and references therein) . Many authors have studied necessary and sufScient conditions 
for the existence of such blow-up time. In particular, in 131 [5l [6l [7] equation 



was considered. In these works, conditions for the existence of a finite blow-up time, as well 
as upper and lower estimates of it, were given, although, in some cases, they were not very 
accurate. 

A way of improving these estimations is to study numerical approximations of the so- 
lution; in this aspect, collocation methods have proven to be a very suitable technique for 
approximating nonlinear integral equations, because of its stability and accuracy (see [5]). 
Hence, the aim of this paper is to check the usefulness of collocation methods for detect- 
ing blow-up solutions of the nonlinear HVHIE (homogeneous Volterra-Hammerstein integral 
equation) given by equation ([T]). 
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In a previous work 0, we used collocation methods for approximating the nontrivial 
solutions of ([T]). In particular we studied the collocation solutions of the implicitly linear 
HVIE (homogeneous Volterra integral equation) 

z{t) = G ((Vz) (t)) ^cl^j^K (t, s) z{s) ds^ , t e /, (2) 

being V the linear Volterra operator, and assuming that some general conditions were held. 
In particular, it was required that the kernel K was locally bounded; nevertheless, non locally 
bounded kernels appear in many cases of blow-up solutions (see [7]) and so, in the present 
work, we shall first replace this condition with another that does not exclude this kind of 
kernels. Hence, the new general conditions that we are going to impose (even if it is not 
explicitly mentioned) are: 

• Over K. The kernel K -.R"^ ^ [0, +oo[ has its support in {(t, s) € : < s < t}. 

For every t > 0, the map s ^ K (t, s) is locally integrable, and IC{t) :— Jq K {t, s) ds is 
a strictly increasing function. 

New condition: limj_^£i+ K {t, s) ds = for all a > 0. 

Note that for convolution kernels (i.e. K{t, s) = k{t ~ s)) this new condition always 
holds, because obviously k is locally integrable and then 

ft f^t — a ft 

lim / k{t — s) ds = lim / k{s) ds ~ lim / k{s) ds = 0. 

t-)-Q+ J a t^a+ Jq t^0+ Jq 

• Over G. The nonlinearity G : [0, +oo[ [0,+oo[ is a continuous, strictly increasing 
function, and G (0) — 0. 

There is a one-to-one correspondence between solutions of (Q^and ([2| (see [HI HO])- Par- 
ticularly, if z is a solution of then y := Vz is a solution of^]). Moreover, y = G^^ o z, 
since G is injective. 

The paper is structured as follows: in Section [2] we introduce briefly the basic definitions 
and state the notation. Next, we devote Section[3]to generalize the results about existence and 
uniqueness of nontrivial collocation solutions obtained in |S| to the new general conditions. In 
Section|4j we introduce the concept of "blow-up collocation solution" and analyse numerically 
some blow-up time estimates using collocation methods in particular examples where the 
previous results about existence and uniqueness can be applied. Finally, in Section [5] we 
discuss the relationships between necessary conditions for blow-up of collocation solutions 
and exact solutions. 



2 Collocation problems for implicitly linear HVIEs 

Following the notation of [8 , a collocation problem for equation ^ is given by a mesh 
Ih ■— {tn ■ — to < ti < . . . < tN — T} and a set of m collocation parameters < ci < 
... < Cm < 1. We denote /i„ := tn+i — tn {n — 0, . . . , N — 1) and the quantity h := 
max{/i„ : < n < — 1} is called the stepsize. Moreover, the collocation points are given 
by in,i := tn + Cihn [i ~ 1, . . . ,m), and the set of collocation points is denoted by Xh (see 
[51[TT]). A collocation solution z^ is then given by the collocation equation 

zh (t) = g(^J^ K {t, s) zh (s) dsj , tGXh, 
where Zh is in the space of piecewise polynomials of degree less than m. 
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Remark 2.1 From now on, a "collocation problem" or a "collocation solution" will be always 
referred to the implicitly linear HVIE So, if we want to obtain an estimation of a solution 
of the nonhnear HVHIE ([I}, then we have to consider yh '■— Vzh- 

Remark 2.2 At the beginning of Section |4| we define the concepts of "blow-up collocation 
problems" and "blow-up collocation solutions" . To do this, we extend the concept of "collo- 
cation solution" to meshes Ih with infinite points. Taking this into account, if zu is a blow-up 
collocation solution, then yh also blows up at the same blow-up time. 

As it is stated in [8] , a collocation solution Zh is completely determined by the coefficients 
Zn.i Zh {tn,i) (n ^ 0, . . . , iV - 1) {i = 1, . . . , to), since Zh (t„ + w/i„) = J2'jLi Lj (v) Zn.j 
for all V e ]0, 1], where Li [v) := 1, Lj {v) := YVk^j c~-Ck (-^ ~ 2, . . . ,to) are the Lagrange 
fundamental polynomials with respect to the collocation parameters. The values of Zji i are 
given by the system 

ZnA = G {tn.i) + /l„ ^ Bn («, j) -^n , (3) 

where ^ 

Bn{i,j)-^ / K {tn,i,tn + Shn) Lj {s) ds 

Jo 

and ^ 

F„ (t) := [" K (t, s) Zh (s) ds. (4) 
Jo 

The term Fn {tn,i) is called the lag term. 

3 Existence and uniqueness of nontrivial collocation so- 
lutions 

Given a kernel K, a nonlinearity G and some collocation parameters {ci, . . . , Cm}, in H] there 
were defined three kinds of existence of nontrivial collocation solutions (of the corresponding 
collocation problem) in an interval / = [0,T] using a mesh //ji existence near zero (col- 
location solutions can always be extended a bit more, but it is not ensured the existence of 
nontrivial collocation solutions for arbitrarily large T), existence for fine meshes (there 
exists a fine enough mesh for which there is a nontrivial collocation solution) , and uncondi- 
tional existence (given any mesh, there always exists a nontrivial collocation solution) . For 
these two last kinds of existence, it is ensured the existence of nontrivial collocation solutions 
in any interval, and hence, there is no blow-up. 

In [5], two particular cases of collocation problems were considered: 

• Case 1: TO = 1 with ci > 0. 

• Case 2: m — 2 with ci = 0. 

In these cases, the system ([3| is reduced to a single nonlinear equation, whose solution is 
given by the fixed points of G (a + /3y) for some a, /3. 

The proofs of the results obtained in ^ rely on the locally boundedness of the kernel 
K. Therefore, we will revise them, imposing the new general conditions and avoiding other 
hypotheses that turn out to be unnecessary. For example, we extend the results to decreasing 
convolution kernels, among others. 
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As in the previous work [S], we wiU need a lemma on nonlinearities G satisfying the 
general conditions: 



Lemma 3.1 The following statements are equivalent to the statement that is unbounded 
(in ]0; +oo[j; 

(i) There exists /3o > such that G {l3y) has nonzero fixed points for all < (5 < /Sq. 

(a) Given A> 0, there exists Pa > such that G {a + j3y) has nonzero fixed points for all 
< a < A and for alio < P < (3a- 

3.1 Case 1: m = 1 with ci > 

First, we shall consider m = 1 with ci > 0. Equations ([S]) are then reduced to 

ZnA = G{Fn{tn,i) + h^B^Z„,i) (n = 0,...,7V-l), (5) 

where ^ 

B„:=S„(1,1)^ / i^(i„,i,t„ + s/i„) ds (6) 
Jo 

and the lag terms Fn [tn.i) are given by Q with i — 1. Note that i?„ > because the 
integrand in Q is strictly positive almost everywhere in ]0,ci[. 

Remark 3.1 In i9i, it is ensured that lim;,^_j.o+ hnBn = taking into account (j6| and the 
old general conditions imposed over K. Thus, we must ensure that the limit also vanishes 
considering the new general conditions. In fact, with a change of variable, we can express ^ 
as ^ 

Bn = ^ [ "'V(i„,i,s) ds. (7) 
Hence, by ^ and the new condition, we have 

lim hnBn = lim / K (t, s) ds ~ 0, 

since = <„ + ci/i„. 

Now, we are in position to give a characterization of the existence near zero of nontrivial 
collocation solutions: 



Proposition 3.1 There is existence near zero if and only if — ^ is unbounded. 



G{v) 

Proof. 

(-4=) Let us prove that if is unbounded, then there is existence near zero. So, we are 
going to prove by induction over n that there exist Hn >0(n = 0,...,A^— 1) such that if 
< /i„ < Hn then there exist solutions of the system (Isj) with Zq.i > 0: 



For n = 0, taking into account Remark 3.1 and Lemma |3. if (^ij (see below), we choose 



a small enough Hq > such that < hoBo < (3o for all < ha < Hq. So, since the lag 
term is 0, we can apply Lemma 3.1-(i) to the equation ([5|, concluding that there exist 
strictly positive solutions for Zq i . 



Let us suppose that, choosing one of those Zq.i and given n > 0, there exist . . . , 
Hn-i > such that if < hi < Hi {i = 1, ... ,n — 1) then there exist coefficients 
Zi 1, . . . , Z„_i 1 fulfilling the equation ([s]). Note that these coefficients are strictly 
positive, and hence, it is guaranteed that the corresponding collocation solution is 
(strictly) positive in ]0,i„]. 
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Finally, we are going to prove that there exists Hn > such that if < /i„ < 
then there exists Z„_i > fulfilling the equation ([5| with the previous coefficients 

^0,1, ■ ■ • Zn^iX. 

Let us define 

A := max{F„ (i„ + cie) : < e < 1} . 

Note that A exists because s ^-^ K {t^ s) is locally integrable for all i > 0, and Zh is 
bounded in [0,t„]; hence, F„ (see Q) is bounded in [i„,i„ +ci]. Moreover, A> Q 
because K and Zh are positive functions. 

So, applying Lemma [3T[ (ii) (see below), there exists > such that G (a + /??/) has 
nonzero (strictly positive) fixed points for a\\ Q < a < A and for all < /3 < /3yi. On 
one hand, taking into account Remark |3.1[ we choose a small enough < < 1 such 
that < hnBn < Pa for all < /i„ < on the other hand, choosing one of those 
hn, we have < F {tn.i) < A. Hence, we obtain the existence of Zn,i as the strictly 
positive fixed point of G {F (tn_i) + hnBny). 



(=>) For proving the other condition, we use Lemma 3.1 -(i), taking into account Remark 



Remark 3.2 A similar result was proved in [5] with the additional hypothesis "K {t, s) < 
K {t', s) for all < s < i < t'" (or "fc is increasing" for the particular case of convolution 
kernels K{t,s) = k{t — s)). We have shown that this hypothesis is not needed and the 
difference between the proofs is the choice of A. 



Taking into account Remark 3.1 we can adapt the results given in [5] about sufficient 
conditions on existence for fine meshes and unconditional existence to our new general con- 
ditions. This will help us to identify collocation problems without blow-up (see Proposition 



3.2 below). But first we need to recall the following definition. 



Definition 3.1 We say that a property V holds near zero if there exists e > such that V 
holds on ]0, 5[ for all < (5 < e. On the other hand, we say that V holds away from zero if 
there exists r > such that V holds on ]t, +00 [ for all t > t. 

The concept "away from zero" was originally defined in a more restrictive way in [S]; 
nevertheless both definitions are appropriate for our purposes. 

Proposition 3.2 Let be unbounded near zero. 

• If K is a convolution kernel and is hounded away from zero, then there is existence 
for fine meshes. 

• // is unbounded away from zero, then there is unconditional existence. 

If, in addition, '^^^ is a strictly decreasing function, then there is at most one nontrivial 
collocation solution. 

Remark 3.3 The proof is analogous to the one given in "9^,. Note that the property 

is unbounded away from zero" appears there as "there exists a sequence {yn}n=i '^^ positive 

real numbers and divergent to +00 such that lim„_>.+oo = 0", but both are equivalent. 
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3.2 Case 2: m = 2 with ci = 

Considering m — 2 with ci = 0, we have to solve the following equations: 

Znj - G(F„(t„a)) (8) 
Zn,2 = G(F„(t„,2) + /i„B„(2,l)Z„,i + /i„B„(2,2)Z„,2), (9) 

for n = 0, . . . , iV — 1, where 

B„(2,j)= / i^(i„,2,t„ + s/i„)Lj (s) ds, (j = l,2) (10) 
Jo 

and F„ (i = 1,2) are given by Q. Note that i?„ (2, j) > for j = 1,2, because the 

integrand in (10) is strictly positive almost everywhere in ]0, C2[. 

Remark 3.4 As in the first case, we have to ensure that lim;i^_j.o+ /i„i?„ (2, j) = for 
j = 1,2, taking into account the new general conditions (see Remark 3.11. Actually, since 
< Lj (s) < 1 for s G ]0, C2[, we have 

0<h„B„{2,j)<K, K{t^^2,tn + shr,) ds = Kitr,a,s) ds. (11) 

Hence, by (11) and the new condition, the following inequalities are fulfilled: 

< lim KBn (2, j) < lim /" K (t, s) ds = 0, 

''"^0+ t^t+ Jt„ 

Analogously to the previous case, we present a characterization of the existence near zero 
of nontrivial collocation solutions: 

Proposition 3.3 Let the map t K {t, s) be continuous in ]s, e[ for some e > and for all 

< s < e (this hypothesis can be removed if C2 ~ 1 ). Then there is existence near zero if and 
only if is unbounded. 

Proof. 

(<^=) Let us prove that if ^^^^ is unbounded, then there is existence near zero. So, we are 
going to prove by induction over n that there exist >0(n = 0,...,A^— 1) such that if 
< /i„ < Hn then there exist solutions of the system (|9|) with Z0.2 > 0: 



• For n = 0, taking into account Remark 3.4 and Lemma 3.1 /ij, we choose a small 
enough Hq > Q such that < /iq^o (2, 2) < /?o for all < /iq < -ffo- So, since the lag 
terms are and Zo,i = G (0) = 0, we can apply Lemma 3.1 /ij to the equation (|9|, 
concluding that there exist strictly positive solutions for ^0,2- 

• Given 71 > 0, let us suppose that, choosing one of those .^0,2, there exist Hi, ... , iJ„_i > 
such that if < hi < Hi {i — 1, ... , n— 1) then there exist coefficients Zi^2, ■ ■ ■ , Zn-1.2 
fulfilling the equation Moreover, let us suppose that Zh is positive in [0,t„], i.e. 
these coefficients satisfy Z; 2 > — C2) Zi i > for / = 1, . . . , n — 1, where Z; 1 is given 
by(§. 

• Finally, we are going to prove that there exists iJ„ > such that if < /i„ < H^ then 
there exists Z„^2 > fulfilling the equation (|9| with the previous coefficients, and Zh is 
positive in [0,t„+i], i.e. Z„_2 > (1 - C2) Z„_i > 0: 

Let us define 

A := max{F„ (t„ + C2e)+ (12) 



eG {F^{tn+cie)) 



K{t^^2,tn + sh„)Li{s) ds : 0<e<l^ 
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Note that A exists and A> because s i-^ K {t, s) is locally integrable for all t > 0, Zh 
is bounded and positive in [0,t„], the nonlinearity G is bounded and positive, and the 
polynomial Li is bounded and positive in [0,C2]. 



So, applying Lemma i.l (ii), there exists > such that G {a + (5y) has nonzero 
(strictly positive) fixed points for all < a < A and for all < /3 < /3yi. On one 
hand, taking into account Remark |3.4[ we choose a small enough < Hn < 1 such 
that < ft.„B„(2,2) < /3a for all < hn < Hn, on the other hand, taking into account 
Q, the lag terms Fn{tn,i) are positive for i = 1,2, because K and Zh are positive. 
Therefore, by ([s]), Zn,i = G(F„(t„.i)) is positive, because G is positive. Moreover, 
hnBn (2, 1) is positive. So, 

< (i„,2) + KBn (2, 1) Zn,i < A. 

Hence, we obtain the existence of Zn.2 as the strictly positive fixed point of the function 
G (a + py) where a := F„ (t„,2) + KB^ (2, 1) Z„,i and f3 := KB„{2, 2). 

Concluding, we have to check that Z„^2 > (1 ^ C2) ^n,i- On one hand, 

Zn,2 = G {t Y,Bn{2,])Zn}^ >G{Fn{tna)) 

because Bn{'2,j)Znj > for j = 1,2. On the other hand, since t 1— >■ K{t,s) is 
continuous in ]s,e[ for some e > and for all < s < e, we can suppose that the 
stepsize is small enough for tn,i < e, and then Fn{t) is continuous in tn,i (see|4]). Hence, 
since G is continuous, G {Fn{t)) is also continuous in i„_i, i.e. lim/j^_^g+ G {Fn (in, 2)) = 
G {Fn {tn,i))- Therefore, choosing a small enough /i„, we have 

G (Fn (i„,2)) > (1 - C2) G (Fn (i„,l)) = (1 - C2) Zn,l. 



(=>) For proving the other condition, we use Lemma 3.1 -(i), taking into account Remark 



Remark 3.5 A similar result was proved in [9 , but the hypothesis on the kernel was 



"if (i, s) < K{t',s) for all < s < t < t'", and so, unlike Proposition 3.3 it could not 
be applied to decreasing convolution kernels (for example). Again, the main difference be- 
tween both proofs lies in the choice of A. 

Note that for convolution kernels K(t,s) — k{t — s), the hypothesis on K is equivalent 
to say that k is continuous near zero, which is a fairly weak hypothesis. On the other hand, 
for general kernels, this hypothesis is only needed for ensuring that the mapping 1 1-^ K (t, s) 
is continuous in tn.i- Therefore, if it is not continuous, it is important to choose a step hi 
such that t i—i' K {t, s) is continuous at the collocation points t^.i (i = 0, . . . — 1). More 
specifically, the hypothesis on K is only needed for checking that Zn,2 > (1 — C2) Zn.i, and 
hence, for ensuring that z/j is positive in [0,t„]. Nevertheless, this hypothesis is not needed 
for verifying that Zo^2 > -^0,1 = 0; and so, it is always ensured the existence of Zi^2 (and 
obviously 2^i,i), even if the hypothesis does not hold. Hence, if we use another method for 
checking that Zi^2 > (1 — C2) Zi^i (for example, numerically), then it is ensured the existence 
of .^2,2- Repeating this reasoning we can guarantee the existence of Zn,2, removing the 
hypothesis on K. 

Following this idea, in the next proposition we state a result analogous to Proposition |3.3| 
without hypothesis on the kernel. 
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Proposition 3.4 is unbounded if and only if there exists Hq > such that there are 

nontrivial collocation solutions in [0,ti] for < Hq < Hq. In this case, there always exists 
Hi > such that there are nontrivial collocation solutions in [0,t2] for < hi < Hi. 

Moreover, if is unbounded and there is a positive nontrivial collocation solution in 
[0,tn], then there exists Hn > such that there are nontrivial collocation solutions in [0,tn+i] 
for 0<hn< Hn. 



Considering Remark 3.4 and taking the coefficient A given in (12 1, we can adapt some 
results of ^ about sufficient conditions on existence for fine meshes and unconditional ex- 
istence to our new general conditions. These results are useful for identifying collocation 
problems without blow-up. 

Proposition 3.5 Let be unbounded near zero, and let the map t i-^ K {t, s) be con- 

tinuous in ]s,e[ for some e > and for all < s < e (this hypothesis can be removed if 
C2 - !)■ 

• If K is a convolution kernel and "^^^ is bounded away from zero, then there is existence 
for fine meshes. 

• // is unbounded away from zero, then there is unconditional existence. 

If, in addition, '^'^ is a strictly decreasing function, then there is at most one nontrivial 
collocation solution. 

Moreover, as in Proposition |3.4[ we can state a result about unconditional existence 
without hypothesis on the kernel: 

Proposition 3.6 Let be unbounded near zero. 

If is unbounded away from zero and there is a positive nontrivial collocation solution 
in \Q,tn], then there is a nontrivial collocation solution in [0,t„+i] for any hn > 0. 

//, in addition, is a strictly decreasing function, then there is at most one nontrivial 

collocation solution. 



3.3 Nondivergent existence and uniqueness 

Our interest is the study of existence of nontrivial collocation solutions using meshes Ih with 
arbitrarily small h > Q. So, we are not interested in collocation problems whose collocation 
solutions "escape" to -l-cxi when a certain hn — > 0+, since this is a divergence symptom. 
Following this criterion, we define the concept of "nondivergent existence" : 

Let = io < . . . < t„ be a mesh such that there exist nontrivial collocation solutions, 
and let Sh„ be the index set of the nontrivial collocation solutions of the corresponding 
collocation problem with mesh to < . . . < t„ < i„ -I- ft,„. Given s € Su^, we denote by Zs-n.i 
the coefficients of the corresponding nontrivial collocation solution verifying equations 
Then, we say that there is nondivergent existence in t^ if 

inf <^ max {Zs.nA \ 

s£Sh„ [_i=l,...,m J 

exists for small enough /i„ > and it does not diverge to +oo when hn — > 0^ . Given a mesh 
Ih = {0 = to < ■ ■ ■ < t^} such that there exist nontrivial collocation solutions, we say that 
there is nondivergent existence if there is nondivergent existence in t^ for n = 0, . . . , N . See 
[S] for a more detailed analysis, where we also define the concept of nondivergent uniqueness. 
In [S| it is proved the next result. 
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Proposition 3.7 In cases 1 and 2 with existence of nontrivial collocation solutions, there is 
nondivergent existence if and only if ^^^^ is unbounded near zero. 

If, in addition, G is "well-behaved" , then there is nondivergent uniqueness. 

We say that G is "well-behaved" if '^^"^^'f is strictly decreasing near zero for all a > 0. 
Note that this condition is very weak (see f9|). 

So, taking into account Propositions |3 . l] |3 . 3| and |3 . 7[ the main result of this paper is: 



Theorem 3.1 (Hypothesis only for case 2 with C2 7^ 1.' the map t 1— > K{t,s) is continuous 
in ]s, e[ for some e > and for all < s < e.) 

There is nondivergent existence near zero if and only if is unbounded near zero. 

If, in addition, G is "well-behaved" then, there is nondivergent uniqueness near zero. 



In the same way, we can combine Propositions |3.4| and |3.7| obtaining a result about 
nondivergent existence and uniqueness without hypothesis on the kernel, and, as Theorem 



3.1 it can be useful to study numerically problems with decreasing convolution kernels in 
case 2 with C2 7^ 1. 

We can also combine Propositions |3. 2[ [33| and [3^ with Proposition [SJj obtaining results 
about nondivergent existence and uniqueness (for fine meshes and unconditional) that are 
useful for identifying collocation problems without blow-up. Some of these results can be 
reformulated as necessary conditions for the existence of blow-up (see Section Isl) . 



4 Blow-up collocation solutions 

In this section we will extend the concept of collocation problem and collocation solution in 
order to consider the case of "blow-up collocation solutions" . 

Definition 4.1 We say that a collocation problem is a blow-up collocation problem (or has 
a blow-up) if the following conditions are held: 

• there exists T > such that there is no collocation solution in / ~ [0, T] for any mesh 
hi, and 

• given M > there exists < r < T, and a collocation solution z^ defined on [0, r] such 
that \zh{t)\ > M for some t e [0,t]. 

We can not speak about "blow-up collocation solutions" in the classic sense, since "collocation 
solutions" are defined in compact intervals and obviously they are bounded; so, we have to 
extend first the concept of "collocation solution" to open intervals / — [0,T[ before we are 
in position to define the notion of "blow-up collocation solution" . 

Definition 4.2 Let / := [0, T[ and Ih be an infinite mesh given by a strictly increasing 
sequence {tn}n=o with to — and convergent to T. 

• A collocation solution on I using the mesh 1^ is a function defined on / such that it is 
a collocation solution (in the classic sense) for any finite submesh {tn}n=o with N €N. 

• A collocation solution on / is a blow-up collocation solution (or has a blow-up) with 
blow-up time T if it is not bounded. 

Given a collocation problem with nondivergent uniqueness near zero, a necessary condition 
for the nondivergent collocation solution to blow-up is that there is neither existence for fine 
meshes nor unconditional existence. So, for example, given a convolution kernel K{t, s) = 
k(t — s), in cases 1 and 2 we must require that is unbounded away from zero (moreover. 
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in case 2 with C2 ^ 1, we must demand that there exists e > such that k is continuous in 

For instance, in |5] it is studied equation (nj) with convolution kernel k{t — s) = {t — s)^ , 
/? > 0, and nonlinearity G{y) = t" , < a < 1, concluding that the nontrivial (exact) solution 
does not have blow-up. If we consider a collocation problem with the above kernel and 
nonlinearity in cases 1 and 2, then we can ensure unconditional nondivergent uniqueness (by 
Propositions 3.2 3.5 and 3.71; thus, we can also conclude that the nondivergent collocation 
solution has no blow-up. Actually, we reach the same conclusion considering any kernel 
satisfying the general conditions (in case 1 and case 2 with C2 = 1) and such that t^K{t, s) 
is continuous in ]s, e[ for some e > and for all < s < e (in case 2 with C2 ^ 1). In Section 
[5] we discuss in more detail the relationships between necessary conditions for blow-up of 
collocation solutions and exact solutions. 



4.1 Numerical algorithms 

First we need to recall the definition of existence near zero, mentioned at the beginning of 
Section [Sj and given in [9]. 

Definition 4.3 We say that there is existence near zero if there exists Hq > such that 
HQ < ho < Hq then there are nontrivial collocation solutions in [0, ti]; moreover, there exists 
Hn > such that if < /i„ < then there are nontrivial collocation solutions in [0, i„+i] 
(for n — 1, . . . , N — 1 and given ho, ... , hn-i > such that there are nontrivial collocation 
solutions in [0, Note that, in general, i?„ depends on ho, . . . , 

Given a blow-up collocation problem with nondivergent uniqueness near zero (always 
in cases 1 and 2), we are going to describe a general algorithm to compute the nontrivial 
collocation solution and estimate the blow-up time. In each iteration, the stepsize hn is 
adapted to Hn, since only if /i„ < Hn we can ensure the existence of the nontrivial collocation 
solution. 

As long as /i„ < Hn, the solution on i„ for each n is obtained from the attracting fixed 
point of a certain function y G{a + I3y), given by equation ^ in case 1, or ([9| in case 2. 
Therefore, the key point in the algorithm is to decide whether there is a fixed point or not 
and, if so, calculate it. 

In order to check that there is no fixed point, the most general technique consists on 
iterate and, if a certain bound (which may depend on the fixed points found in the previous 
steps) is overcomed, then it is assumed that there is no fixed point. 

So, if for a certain n there is no fixed point, then a smaller hn should be taken (e.g. hn/2 
in the examples below), and this procedure is repeated with the new a, /3 corresponding to 
the new h„. If /i„ becomes smaller than a given tolerance (e.g. 10""'^^ in our examples), the 
algorithm stops and — ho + . . . + hn is the estimation of the blow-up time. 

Nevertheless, in particular cases, specific algorithms can also be used, as in the case of 
Example 1, where the fixed points can be found analytically: 

• If < a < 1: 

- If 1 - a < /? < ^ the fixed point is 

- If /3 < 1 - a then the fixed point is ^ (^l3 + ^ fi"^ 4- 4a^ . 

• If a > 1 and ^ < ^, then the fixed point is 

If /3 > ^ then there is no solution for the corresponding hn- Hence a smaller /i„ is taken 
as described in the general case. 

However, both the computational cost and accuracy of the general and specific algorithms 
are quite similar. 
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4.2 Examples 



We will consider the following examples 
1. K{t,s) = 1 and G{y) 



Vy if ye [0,1] 

ifye]l,+oo[ 



2. K {t, s) = 1 and G (y) 



if 2/ e [0, 1] 
if y G ]1, +oo[ 



3. K (t, s) = k{t - s) = t - s and G (y) of Example 1. 

4. K {t, s) = k{t - s) = {-K{t- s)y^'^ and G (y) of Example 1. 

In the first two examples, the blow-up time of the corresponding (exact) solution is i = 3, 
while in the other two the blow-up time is unknown. In [7 it was studied a family of equations 
to which Example 4 belongs. 

In all the examples, both K and G fulfill the general conditions and the hypotheses of 
Theorem |3.1[ and thus it is ensured the nondivergent existence and uniqueness near zero 
(note that since the kernel in Example 4 is decreasing and unbounded near zero, this case 
can not be studied in [5]). Moreover, G(y)/y is unbounded away from zero and y/G(y) is 



bounded away from zero, and hence there can exist a blow-up (see Propositions 3.2 and 3.5 ). 

He have found the numerical nondivergent collocation solutions for a given stepsize h 
(= /iq), using the algorithms described in Sectio n |4.1| (the specific for the Example 1 and 
the general for the rest). In Figures [T] [2j [3] and |4| the estimations of the blow-up times of 
the collocation solutions for a given stepsize h are depicted, varying ci in case 1 or C2 in 
case 2. Note that the different graphs for different h intersect each other in a fairly good 
approximation of the blow-up time, as it is shown in Tables [l] [2] |4] and [5| 

In Examples 3 and 4 we do not know the exact value of the blow-up time. However, 
in order to make a study of the relative error analogous to the previous examples, we have 
taken as blow-up time for Example 3 the approximation for h = 0.001 in case 2 with C2 = 0.5: 
t = 5.78482; and for Example 4 the approximation for h = 0.0001 in case 2 with C2 = 2/3 



[Radau I collocation points): t = 1.645842. Results are shown in Tables [7j|8 10 and 11 

The relative error varying ci (case 1) or C2 (case 2) is the "relative vertical size" of the 
graph, and it decreases at the same rate as h. On the other hand, the relative error of the 
intersection decreases faster, sometimes at the same rate as . Moreover, in case 1, the best 
approximations are reached with ci « 0.5, and in case 2 with C2 ~ 2/3 (approximately Radau 



I) for Examples 1, 2 and 4 (see Tables |3j [6] and 12), while for Example 3 the best approxi- 
mations are reached with C2 ~ 0.5; however, in Table |9] are also shown the approximations 
and their corresponding errors for the Radau I collocation points. On the other hand, the 
intersections method offers better results, but at a greater computational cost. 



5 Discussion and comments 

The main necessary conditions for a collocation problem to have a blow-up are obtained from 



Propositions 3.2 (case 1) and 3.5 (case 2), and are mostly related to the nonlinearity, since 
the assumption on the kernel "the map 1 1-^ K (t, s) is continuous in ]s, e[ for some e > and 
for all < s < e" is only required in case 2 with C2 1 and it is a very weak hypothesis. 
Hence, assuming that the kernel satisfies this hypothesis, the main necessary condition for 
the existence of a blow-up is: 

1. is bounded away from zero. 
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In addition, for convolution kernels, there is another necessary condition: 

2. is unbounded away from zero, i.e. there exists a sequence {yn}n^i positive real 

numbers and divergent to +00 such that lim„_j.+oo ofy ) — 0. 

In [4] it is given a necessary and sufficient condition for the existence of blow-up (exact) 
solutions for equation ([l]) with a kernel of the form K(t, s) = (t — s)"^^r(s) with a > 0, r 
nondecreasing and continuous for s 7^ 0, r(s) = for s < 0, and r(s) > for s > 0: 

(— ) ^^<+„o, ^>0^ (13) 

This also holds for convolution kernels of Abel type K{t, s) — {t — s)"~^ with a > 0, gener- 
alizing some results given in [3 . Next we will show that necessary conditions 1 and 2 above 



mentioned are also necessary conditions for the integral given in (13) to be convergent, and 
thus for the existence of a blow-up (exact) solution. 



Proposition 5.1 // (13) holds, then is unbounded away from zero 



Proof. Let us suppose that ^^^^ is bounded away from zero. So, there exists M > such 
that ^ <M for all y > S. Hence 

+00/ „ \ ds ^ /^y/" /■+°° ds _ 
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G{s)J s \M J Js s 



Proposition 5.2 // {IS) holds, then -q^^ is bounded away ft 



rom zero. 



Proof. By Proposition |5.l| is unbounded away from zero and hence, there exists a 

strictly increasing sequence {yn}n=i with yi > S and divergent to -l-oo such that ^ < ^ 
for all n. 

Let us suppose that is unbounded away from zero; so, we can choose {yn}n=i such 
that there exist j/Jj e ]yn,yn+i[ with tjt^ = 1 for each n. Moreover, since the nonhnearity 



G is positive and strictly increasing, we have that ^ > ^ = 1, and then y'^ > G (j/n). 
Hence, we have 



l/a\ / / \ l/a\ 

yn\ \ ^ / -, / yn \ \ ^ a 



Therefore, ( 13 1 does not hold. 



An interesting problem which is not fully resolved is to continue clarifying the relationship 
between the existence of blow-up in exact solutions and in collocation solutions. 



12 



References 



[1] C. M. Kirk. Numerical and asymtotic analysis of a localized heat source undergoing 
periodic motion. Nonlinear Anal. 71 (2009), e2168-e2172. 

[2] F. Calabro, G. Capobianco, Blowing up behavior for a class of nonlinear VIEs connected 
with parabolic PDEs. J. Comput. Appl. Math. 228 (2009), 580-588. 

[3] W. Mydlarczyk, A condition for finite blow-up time for a Volterra integral equation. J. 
Math. Anal. Appl. 181 (1994), no. 1, 248-253. 

[4] W. Mydlarczyk, The blow-up solutions of integral equations. Colloq. Math. 79 (1999), 
no. 1, 147-156. 

[5] T. Malolcpszy, W. Okrasihski, Conditions for blow-up of solutions of some nonlinear 
Volterra integral equations. J. Comput. Appl. Math. 205 (2007), 744-750. 

[6] T. Malolepszy, W. Okrasinski, Blow-up conditions for nonlinear Volterra integral equa- 
tions with power nonlinearity. Appl. Math. Letters. 21 (2008), 307-312. 

[7] T. Malolepszy, W. Okrasihski, Blow-up time for solutions for some nonlinear Volterra 
integral equations. J. Math. Anal. Appl. 366 (2010), no. 1, 372-384. 

[8] H. Brunner, Collocation Methods for Volterra Integral and Related Functional Differen- 
tial Equations, Cambridge University Press, Cambridge, 2004. 

[9] R. Benitcz, V. J. Bolos, Existence and uniqueness of nontrivial collocation solutions of 
implicitly linear homogeneous Volterra integral equations. J. Comput. Appl. Math. 235 
(2011), no. 12, 3661-3672. 

[10] M. A. Krasnosel'skii, P. P. Zabreiko, Geometric Methods of Nonlinear Analysis, Springer 
Verlag, New York, 1984. 

[11] H. Brimncr, Implicitly linear collocation methods for nonlinear Volterra integral equa- 
tions, Appl. Numer. Math. 9 (1992), 235-247. 



13 



A Figures and tables 




Figure 1: Example 1. Numerical estimation of the blow-up time of collocation solutions vary- 
ing ci (case 1) or C2 (case 2), for different stepsizes h. The blow-up time of the corresponding 
(exact) solution is f = 3. 
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Case 1 

(m = 1, ci > 0) 


Varying ci 


Intersection 


min. t 


max. t 


Rel. error 


t 


Rel. error 


h = 0.1 


2.60 

(ci = 1) 


3.66 
(ci = 0.01) 


4- 10-1 


3.03 
(ci = 0.39) 


10-2 


h = 0.05 


2.78 

(Cl = 1) 


3.40 
(ci = 0.01) 


2 • 10-1 


h = 0.01 


2.92 

(ci = 1) 


3.11 
(ci = 0.01) 


6-10-2 


3.002 
(ci = 0.44) 


7 - 10-4 


h = 0.005 


2.96 

(ci = 1) 


3.06 
(ci = 0.01) 


3-10-2 



Table 1: Example 1. Numerical data of Figure [T] (case 1). 



Case 2 

(to = 2, ci 0) 


Varying C2 


Intersection 


min. t 


max. t 


Rel. error 


t 


Rel. error 


h = 0.1 


2.95 

(C2 = 1) 


3.20 

(C2 = 0.01) 


8 - 10-2 


3.0003 

(c2 = 0.6178) 


10-4 


h = 0.05 


2.98 

(C2 = 1) 


3.10 

(C2 = 0.01) 


4- 10-2 


/i = 0.01 


2.995 

(C2 = 1) 


3.020 

(C2 = 0.01) 


8 - 10-3 


3.000003 
(c2 = 0.6232) 


10-6 


h = 0.005 


2.998 

(C2 = 1) 


3.010 

(C2 = 0.01) 


4 - 10-3 



Table 2: Example 1. Numerical data of Figure [T] (case 2). 





Case 1 - ci — 0.5 


Case 2 - Radau I 


Blow-up 


Rel. error 


Blow-up 


Rel. error 


/i = 0.1 


2.933883 


2.2 - 10-2 


2.995253 


1.6 - 10-3 


h = 0.05 


2.965002 


1.2 - 10-2 


2.997602 


8 - 10-4 


= 0.01 


2.992885 


2Aq ■ 10-3 


2.999519 


1.6 - 10-4 


h = 0.005 


2.996434 


1.2 - 10-3 


2.999759 


8 - 10-'"^ 



Table 3: Example 1. Numerical estimations and relative errors of the blow-up time of 
collocation solutions in case 1 with ci = 0.5, and case 2 with Radau I collocation points, 
ci = 0, C2 = 2/3, for different stepsizes h. 
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Casel 



Case 1 




Figure 2: Example 2. Numerical estimation of the blow-up time of collocation solutions vary- 
ing ci (case 1) or C2 (case 2), for different stepsizes h. The blow-up time of the corresponding 
(exact) solution is t = 3. 
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Case 1 

(to = 1, ci > 0) 


Varying ci 


Intersection 


min. t 


max. t 


Rel. error 


t 


Rcl. error 


h = 0.1 


2.68 

(ci = 1) 


3.50 
(ci = 0.01) 


3 • 10-1 


3.015 
(ci = 0.399) 


5 • 10-3 


h = 0.05 


2.82 

(ci = 1) 


3.30 
(ci = 0.01) 


1.5 • 10-1 


h = 0.01 


2.94 

(ci = 1) 


3.08 
(ci = 0.01) 


4.6 • 10-2 


3.0008 
(ci = 0.444) 


2.6- 10-* 


h = 0.005 


2.98 

(Cl = 1) 


3.04 
(ci = 0.01) 


2.3 • 10-2 



Table 4: Example 2. Numerical data of Figure [2] (case 1). 



Case 2 

(to = 2, Ci = 0) 


Varying C2 


Intersection 


min. t 


max. t 


Rel. error 


t 


Rel. error 


h = 0.1 


2.97 

(C2 = 1) 


3.10 

(C2 = 0.01) 


4.3 • 10-2 


3.00012 
(c2 = 0.6466) 


4- 10-5 


0.05 


2.99 

(C2 = 1) 


3.05 

(C2 = 0.01) 


2 • 10-2 


h = 0.01 


2.997 

(C2 = 1) 


3.010 

(C2 = 0.01) 


4.3 • 10-3 


3.000011 
(c2 = 0.65114) 


3.6-10-'^ 


= 0.005 


2.9985 

(C2 = 1) 


3.005 

(C2 = 0.01) 


2.15 • 10-3 



Table 5: Example 2. Numerical data of Figure [2] (case 2). 





Case I - ci — 0.5 


Case 2 - Radau I 


Blow-up 


Rel. error 


Blow-up 


Rel. error 


h = 0.1 


2.941795 


2 • 10-2 


2.999559 


1.5 • 10-* 


h = 0.05 


2.970341 


10-2 


2.999778 


7.4 • 10-5 


/i = 0.01 


2.993979 


2 • 10-3 


2.999955 


1.5 • 10-5 


h = 0.005 


2.996983 


10-3 


2.999978 


7.3 • 10-6 



Table 6: Example 2. Numerical estimations and relative errors of the blow-up time of 
collocation solutions in case 1 with ci = 0.5, and case 2 with Radau I collocation points, 
ci = 0, C2 = 2/3, for different stepsizes h. 
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Casel 



Casel 
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Case 1 

(to = 1, Ci > 0) 


Varying ci 


Intersection 


niin. t 


max. t 


Rcl. error 


t 


Rel. error 


h = 0.1 


5.20 

(ci = 1) 


6.55 
(ci = 0.01) 


2.3- 10-1 


5.8055 
(ci = 0.3961) 


3.6-10-3 


h = 0.05 


5.45 

(ci = 1) 


6.20 
(ci = 0.01) 


1.3-10-1 


h = 0.01 


5.71 

(ci = 1) 


5.89 
(ci = 0.01) 


3.1 - 10-2 


5.7856 
(ci = 0.441) 


1.4- 10-4 


h = 0.005 


5.735 

(ci = 1) 


5.840 
(ci = 0.01) 


1.8- 10-2 



Table 7: Example 3. Numerical data of Figure [3] (case 1). The blow-up time of the corre- 
sponding (exact) solution is assumed to be i « 5.78482. 



Case 2 

(to = 2, Ci = 0) 


Varying C2 


Intersection 


min. t 


max. t 


Rel. error 


t 


Rel. error 


h 0.1 


5.67 

{C2 = 1) 


6.10 

(c2 = 0.01) 


7.4- 10-2 


5.78304 
(c2 = 0.5286) 


3 - 10-4 


h = 0.05 


5.73 

(C2 = 1) 


5.95 
(C2 = 0.01) 


3.8-10-2 


h = 0.01 


5.775 

(C2 = 1) 


5.810 
(c2 = 0.01) 


6 - 10-3 


5.7848 
(C2 = 0.514) 


3.5 - 10-6 


h = 0.005 


5.779 

(C2 = 1) 


5.800 
(c2 = 0.01) 


3.6-10-3 



Table 8: Example 3. Numerical data of Figure [3] (case 2). The blow-up time of the corre- 
sponding (exact) solution is assumed to be i w 5.78482. 





Case 1 - ci = 0.5 


Case 2 - Radau I 


Case 2 - C2 = 0.5 


Blow-up 


Rel. error 


Blow-up 


Rel. error 


Blow-up 


Rel. error 


h = 0.1 


5.694865 


1.6-10-2 


5.751797 


5.7-10-3 


5.788215 


5.9-10-4 


h = 0.05 


5.739117 


7.9-10-3 


5.767963 


2.9- 10-3 


5.786612 


3.1 - 10-4 


h = 0.01 


5.775121 


1.7-10-3 


5.781037 


6.5 - 10-4 


5.784995 


3 - 10-=^ 


h = 0.005 


5.780132 


8.1 - 10-4 


5.783066 


3 - 10-4 


5.784875 


9.5-10-6 



Table 9: Example 3. Numerical estimations and relative errors of the blow-up time of 
collocation solutions in case 1 with ci = 0.5, and case 2 with Radau I collocation points, 
Ci = 0, C2 = 2/3, and ci = 0, C2 = 0.5, for different stepsizes h. The blow-up time of the 
corresponding (exact) solution is assumed to be t w 5.78482. 
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Casel 



Casel 




20 



Case 1 

(m = 1, ci > 0) 


Varying ci 


Intersection 


mill, t 


max. t 


Rel. error 


t 


Rel. error 


h = 0.01 


1.593 

(ci = 1) 


1.722 
(ci = 0.05) 


7.8 • 10-2 


1.6467 

(ci = 0.444) 


5.2 • 10-4 


h = 0.005 


1.613 

(ci = 1) 


1.691 
(ci = 0.05) 


4.7- 10-2 


h = 0.001 


1.6367 

(ci = 1) 


1.6571 
(ci = 0.05) 


1.2 • 10-2 


1.64595 
(ci = 0.461) 


6.6- 10-5 


h = 0.0005 


1.6406 

(ci = 1) 


1.6521 
(ci = 0.05) 


7 • 10-3 



Table 10: Example 4. Numerical data of Figme|4] (case 1). The blow-up time of the corre- 
sponding (exact) solution is assumed to be t « 1.645842. 



Case 2 

(m = 2, ci = 0) 


Varying C2 


min. i 


max. t 


Rel. error 


h = 0.01 


1.6445 
(c2 = 0.95) 


1.6490 
(c2 = 0.33) 


2.7 • 10-3 


h = 0.005 


1.6456 
(c2 = 0.9) 


1.6476 
(c2 = 0.3) 


1.2 • 10-3 


h = 0.001 


1.64577 
(c2 = 0.9) 


1.64622 
(c2 = 0.3) 


2.7-10-4 


h = 0.0005 


1.64581 

(c2 = 0.85) 


1.64607 
(C2 = 0.28) 


1.6 • 10-4 



Table 11: Example 4. Numerical data of Figure |4] (case 1). The blow-up time of the corre- 
sponding (exact) solution is assumed to be t « 1.645842. 





Case 1 - ci = 0.5 


Case 2 - Radau I 


Blow-up 


Rel. error 


Blow-up 


Rel. error 


ft = 0.01 


1.638700 


4.3- 10-3 


1.646991 


7 • 10-4 


ft = 0.005 


1.642843 


1.8- 10-3 


1.646540 


4.2- 10-4 


ft = 0.001 


1.645172 


4.1 • 10-4 


1.645985 


8.7-10-5 


ft = 0.0005 


1.645491 


2.1 • 10-4 


1.645919 


4.7- 10-5 



Table 12: Example 4. Numerical estimations and relative errors of the blow-up time of 
collocation solutions in case 1 with ci = 0.5, and case 2 with Radau I collocation points, 
ci — 0, C2 — 2/3, for different stepsizes ft. The blow-up time of the corresponding (exact) 
solution is assumed to be t « 1.645842. 
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